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A LUMINOUS RECURRENT SUPERSOFT X-RAY SOURCE IN NGC 300 
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ABSTRACT 

We report the results of XMM-Newton observations for an especially luminous supersoft X-ray source (SSS) 
with bolometric luminosity of 10^^ erg s~' in the spiral galaxy NGC 300. The source was detected as a SSS in 
1992 and disappeared in subsequent X-ray observations. The source was active again during recent XMM-Newton 



observations. It appeared to be very soft (kT ~ 60 eV) and very luminous {^^ 10 -10 erg s" ). The two XMM- 
Newton observations also reveal that the source went from a "high" state to a "low" state in 6 days. We also 
found a 5.4-hr periodicity during the "low" state. We consider white dwarf, black hole, and neutron star models 
to explain the nature of the source. 

Subject headings: galaxies: individual (NGC 300) — X-rays: binaiies — X-rays: galaxy 



1. INTRODUCTION 

Supersoft X-ray sources (SSSs) form a distinct class of ob- 
jects, first established through ROSAT observations. The hall- 
marks of SSSs are very soft X-ray emission (kT typically < 100 
eV) and bolometric luminosities of lO-'^- lO-'** erg s"'. The ad- 
vent of Chandra and XMM-Newton provides good opportuni- 
ties to detect and study SSSs in nearby galaxies. Luminous 
|-jQ38-4{)gj.g g-i-) have been found in several nearby galax- 
ies including M31 (Kong et al. 2002; DiStefano et al. 2002a), 
M81 (Swai-tz et al. 2002), MlOl (Pence et al. 2001; Mukai et al. 
2003; DiStefano & Kong 2003 [DK03]), NGC4697 (Sarazin, 
Irwin & Bregman 2001; DK03), M51 (DK03), and M83 (Soria 
& Wu 2003; DK03). 

NGC 300 is an SA(s)d galaxy, seen near face-on (inclination 
angle 46°; Tully 1988) at a distance of 2.0 ± 0.1 Mpc (Freed- 
man et al. 2001). The galaxy has been observed by ROSAT 
(Read et al. 1997; Read & Pietsch 2001) and XMM-Newton 
(Soria & Kong 2003). 

In this Letter, we report the reappearance of a luminous SSS 
in NGC 300 as observed with XMM-Newton. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. ROSAT 

NGC 300 was observed by ROSAT five times from 1991 to 
1997. A detailed analysis of the ROSAT data was done by Read 
& Pietsch (2001). Briefly, the datasets consist of 2 Position 
Sensitive Proportional Counter (PSPC) and 3 High Resolution 
Imager (HRJ) pointings. The exposures range from ~ 9 ks to 
~ 37 ks. The luminous supersoft source was only detected in a 
37 ks PSPC observation taken on 1992 May and June (see Read 
& Pietsch 2001). Spectral analysis was also done by Read et al. 
(1997). We have reanalyzed the PSPC spectrum taken on 1992 
May/June, and used other data to set upper limits on the lumi- 
nosity for the long-term lightcurve. 

We extracted the source spectrum from a 30" circular region, 
while background was from an annulus region (45" and 60" 
radii) centered on the source. The spectrum was grouped into 
at least 20 counts per spectral bin to allow statistics to be 
used. 

2.2. XMM-Newton 



The XMM-Newton instrument modes were full-frame, 
medium filter for the three European Photon Imaging Cameras 
(EPIC). The first observation was taken on 2000 December 26 
for 37 ks and the second observation was on 2001 January 1 
for about 47 ks. After rejecting intervals with a high back- 
ground level, we considered a good time interval of ^ 28 ks 
and ^ 40 ks for the first and second observation, respectively. 
Data were reduced and analyzed with the XMM-Newton SAS 
package v5.3.3. 

We used here the MOS images to determine the position 
of the source because the spatial resolution of MOS detec- 
tor is slightly better than that of pn. The source is located at 
R.A.=00h55ml0s.7,Dec.=-37°38'55" (J2000), -4' (2.4 kpc) 
from the galactic center; the derived positions from the two 
XMM-Newton observations agree with each other This position 
is about 5" off from previous ROSAT PSPC observation, and is 
consistent with the positional error (7.3") quoted by Read & 
Pietsch (2001). 

Source spectra and lightcurvesof XMMU J005510.7-373855 
were extracted with the SAS task XMMSELECT. Source-free re- 
gions were used for background to avoid the chip boundary and 
a nearby faint source. In order to allow statistics to be used, 
all the spectra were binned such that there are at least 20 counts 
per spectral bin. Response matrices were created by RMFGEN 
and ARFGEN. 



3. ANALYSES AND RESULTS 

Spectral analysis was performed by making use of XSPEC 
vl 1 .2. Table 1 shows the best-fitting spectral parameters for the 
three ROSAT and XMM-Newton observations. 

For the ROSAT observation, the spectrum (see Figure 1) 
can be fit with a blackbody model with Nh = 1.1 x 10^' cm"^ 
and kT = 48.7 eV; the 0.2-2 keV luminosity is 10^^'' ergs"'. 
We note that Read et al. (1997) fit the spectrum with a 
thermal bremsstrahlung model with temperature of a 0. 1 keV 
and we confirmed that it is also an acceptable model. Us- 
ing the thermal bremsstrahlung model, the luminosity becomes 
3.8 X lO-'^ ergs"' [the luminosity quoted by Read et al. (1997) 
is 1.6 X lO-'^ erg s"', which is an absorbed luminosity corrected 
for absorption in our own Galaxy], a factor of 2.5 lower than 
the blackbody model. While we cannot distinguish between 
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the two models statistically, we prefer the blackbody model 
as subsequent XMM-Newton observations confirm the super- 
soft nature of the source (see below). Also, the XMM-Newton 
observations would have detected the high energy photons as- 
sociated with the thermal bremsstrahlung model, had the model 
been correct. Finally, the XMM-Newton observations provide 
more photons than ROSAT. 

We fit the pn, MOSl and MOS2 data simultaneously with 
several single-component models with interstellar absorption 
(including absorbed power-law, thermal bremsstrahlung, black- 
body and Raymond- Smith); only the blackbody model provides 
an acceptable fit. The blackbody temperature of both observa- 
tions ranges between 57 eV and 67 eV, while the Nh varies from 
1 .5 X 10^' cm"^ in the first observation to 5.3 x 10^" cm~^ in the 
second observation. The 0.2-2 keV luminosity also drops from 
9 X 10^^ ergs"^ to 1.2 x 10^^ ergs"\ indicating that the source 
is a variable on timescale of days. The spectra of the two ob- 
servations are shown in Figure 1 . 

The combined ROSAT and XMM-Newton long-term 
lightcurve of the SSS is shown in Figure 2, which is con- 
structed from a series of ROSAT and XMM-Newton pointings. 
The source was below the detection limit of other ROSAT ob- 
servations. In these cases, we estimated the 3(t limits, assuming 
a blackbody model with mean Nh (10^^ cm"^) and kT (58 eV). 
The value of Nh plays an important role in estimating the lu- 
minosity. For instance, if we lower the Nh to 5.3 x 10^" cm"^, 
as found in the second XMM-Newton observation, the upper 
limits decrease by a factor of ^ 3. In other words, the source 
varies by as much as factor of 30 between the "low" state and 
the "high" state spanning 8.5 years. 

To search for short-term variability, we computed the Lomb- 
Scargle periodogram (LSP; Lomb 1976; Scargle 1982), a mod- 
ification of the discrete Fourier transform which is generalized 
to the case of uneven spacing. In each observation, we extracted 
the combined background subtracted Ughtcurve from pn and 
MOSs to increase the signal-to-noise. Individual lightcurves 
from each of the three detectors were also used to verify the 
result. By applying the LSP to the combined lightcurve, we 
found that there is a sharp peak at 5.4 hr in the "low" state ob- 
servation (see Figure 3). Independent checks from the pn and 
MOS data also confirmed the periodicity. We determined the 
99.9% confidence level by generating Gausiaan noise datasets 
with the same time intervals and variance as the true data, and 
then performed the LSP analysis on the resulting datasets (see 
Kong, Charles & Kuulkers 1998). The peak power in each 
periodogram (which must be purely due to noise) was then 
recorded. This was repeated 10000 times to obtain good statis- 
tics. The peak at 5.4 hr is well above the 99.9% confidence 
level. The folded lightcurve of the "low" state data in 5.4 hr is 
also shown in Figure 3. 

Similar analysis was also carried out in the "high" state data, 
but there is no significant peak in the LSP. The folded lightcurve 
shows no obvious periodic variabiUty at 5.4 hr (see Figure 3). 
It is not clear if it is due to geometric effect but the exposure 
time of the "high" state covers less than 8 hr, corresponding to 
1.4 cycles of the 5.4-hr period in contrast to the 2.1 cycles in 
the "low" state. We therefore performed Kolmogorov-Smimov 
test to test for variability and found significant variability at the 
99.9% levels. We also search for variability of the hardness 
ratio but no significant change was found in both observations. 

We also examined an Digitized Sky Survey image of NGC 

^ see http://www.aip.de/~jcg/sss/ssscat.html for updated catalog 



300 to search for possible optical counterpart. Within a 4" er- 
ror circle of the XMM-Newton position, there is no obvious op- 
tical counterpart with limiting B and R magnitudes of 21 and 
18.1, respectively. However, the source is surrounded by Hii 
regions (Deharveng et al. 1988; Blair & Long 1997), OB asso- 
ciations (Pietrzynski et al. 2001) and supernova remnants (Blair 
& Long 1997). The nearest Hii region (source 136; Deharveng 
et al. 1988) is located about 20" (~ 200 pc) south-east of the 
X-ray source. Within 50" (~ 500 pc) of the SSS, we found 
an OB association (AS99; Pietrzynski et al. 2001), two super- 
nova remnant candidates (N300-S21, N300-N23; Blair & Long 
1997) and a Hii region (N300-H20; Blair & Long 1997). We 
also looked for UV counterpart from the Optical Monitor (OM) 
image (see Soria & Kong 2003 for details of the OM observa- 
tions); no counterpart is found at the position of the SSS in both 
observations. 

4. DISCUSSION 

It is clear from the spectral fits of all ROSAT and XMM- 
Newton observations that the XMMU J0055 10.7-373855 is 
very soft {kT ~ 48 - 67 eV) and is highly variable. Such a soft 
spectrum is consistent with SSSs seen in our own Galaxy, the 
Magellanic Clouds, M31 (Greiner 2000 ^) and several nearby 
galaxies (see § 1). Since the Galactic and Magellanic Cloud 
SSSs are conjectured to be white dwarfs (WDs), it is reason- 
able to consider a WD model for the source in NGC 300. Since, 
however, the luminosity in the high state may be too large to be 
consistent with the standard SSS WD model, we also consider 
black hole (BH) and neutron star (NS) models. 

The X-ray emission of the SSS could come from a WD that 
burns accreting hydrogen in a quasi-steady manner. The max- 
imum luminosity is the Eddington luminosity for a 1.4M0 ob- 
ject. WD models therefore seem to be ruled out by the high 
bolometric luminosities during 2 of the observations. It is nev- 
ertheless worth noting that the WD model temperatures are 
comparable to those we observed, and that X-ray variabiUty 
by a factor of a few over times of days and months has been 
observed in SSSs. (see e.g., Greiner and Di Stefano 2002 and 
references therein.) If the 5.4-hour variation detected during the 
X-ray low state corresponds to an orbital period, then this sys- 
tem may be similar to 1E0035. 4-7230, a Magellanic Cloud SSS 
that appears to have a 4 hour period (Kahabka 1996). In such 
a short-period binary, it is likely that a radiation-driven wind 
plays an important role in the binary evolution, (see, e.g., van 
Teeseling & King 1998) A WD model would perhaps be the 
most natural explanation for the SSS in NGC 300, were it not 
for the high estimated bolometric luminosities. Note that it may 
be difficult to use beaming to circumvent the Eddington limit in 
these systems, since the inner region of the accretion disk can 
also become a luminous emitter of SSS radiation (Popham & 
Di Stefano 1996). 

SSS emission is expected from accreting BHs. Modeling 
the accretion disk as a thin disk which is optically thick, we 
can derive a relationship between the minimum mass of the 
accretor, and the observed temperature and luminosity (Frank, 
King, & Raine 2002). Using kT = 60 eV and Lx = 10^^ erg s'' 
(10^^ erg s"'), we find that the accretor mass is greater than 
« 890Mq (2800Mq). Although the fact that the observed lu- 
minosity is lO'^ — 10^ times smaller than the Eddington lumi- 
nosity may call the validity of the disk model into question, this 
calculation suggests that the accretor is an intermediate-mass 
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BH. If the accretor is a BH, we may use Lots = Q.\mc^{rj/Q.\), 
where m is the accretion rate and 77, the efficiency factor, is 
likely to be close to 0.1. This yields m « 1.8 x lO"**M0/yr 
(1.8 X IQ-'^Mo/yr) for Uh, = lO^'^ergs"' (Lohs = lO^^ergs"'). 
This would be consistent with an irradiation-driven wind from a 
low-mass donor. Measurement of the orbital period could pro- 
vide supporting evidence. Indeed, if the donor fills its Roche 
lobe, then Pgrb ~ (8.9 hrs)iMd/MQ), where is the mass of 
the donor star. If P„rh = 5.4 hrs, then « 0.61Mq. 

It is of course possible that the observed emission emanates 
from an accreting NS. In this case the SSS emission would pre- 
sumably emanate from a photosphere which is much larger than 
the neutron star itself. The photospheric radius would be differ- 
ent in each observation; 3.7 x 10'' cm during the ROSAT oh- 
servation, 2.6 x 10^ cm during the XMM-Newton high-state ob- 
servation, and 6.8 x 10^ cm during the XMM-Newton low-state 
observation, corresponding to the different temperatures and lu- 
minosities. 

At this point, neutron star models seem conceptually 
unattractive because there is no obvious explanation for why 
the photosphere should achieve these relatively large sizes, and 
no way to relate them to the system's fundamental physical pa- 
rameters. In addition, just as in the WD models, the luminos- 
ity appears to be super-Eddington during two of the observa- 
tions, unless the neutron star has a mass as large as Z-AMq. 
In this case, if the efficiency of turning accretion energy into 
X-rays is the same during all 3 observations, then the accretion 



rate must have changed by an order of magnitude, approaching 
2 X 1O"**M0 yr"' during the high state, (see DiStefano et al. 
2002b for a similar model of the X-ray source in the globular 
cluster Bo 375.) If, on the other hand, the mass of the neutron 
star is close to \AMq, then other ways around Ihc Eddington 
limit must be found; e.g., the energy could be beamed into a 
smaller solid angle during the 2 high-state observations. In this 
case, the accretion rate could be significantly lower. 

5. SUMMARY 

We have studied a recurrent luminous SSS in NGC300 with 
XMM-Newton. The source was seen by ROSAT in 1992 and 
fell below the detection Umit in subsequent ROSAT observa- 
tions. It reappeared in recent XMM-Newton observations with 
bolometric luminosities between 10-"^ and lO-''' ergs"'. During 
the "low" state, the source showed a 5.4-hr periodicity. If we 
consider the periodicity is due to the orbital period, then the 
X-ray emission can be explained by WD, BH, and NS models. 
However, WD and NS models appear to be unlikely due to the 
high X-ray luminosity during the "high" state. Eurther repeated 
X-ray observations of the SSS in different states may permit 
discrimination among these models. 

This work was supported by NASA under an LTSA grant, 
NAG5-10705. A.K.H.K. acknowledges support from the 
Croucher Foundation. 
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Table 1 

Best-htting Spectral Parameters 



Date 


Nh 

(xlO" cm"2) 


kT 
(eV) 


Lx 


Lboi 


" x'/d.o.f 


1992 May 26 
2000 Dec 26 
2001 Jan 1 


i.08:!SJi 

0.53!g;§? 


48.7!^'l 
56.8li 
66.5:i;^ 


10 

9 

1.1 


20 
15 
1.9 


0.74/11 
1.07/39 
1.15/48 



Note. — All quoted uncertainties are Icr. 

^0.2-2 keV luminosity ( x lO'* erg s"'), assuming a distance of 2 Mpc 
''Bolometric luminosity (x 10^* erg s~') 



ROSAT PSPC (1992 May/June) XMM-Newton (2000 December) 
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Fig. 1.— Energy spectra of XMMUJ005510.7-373855 (see Table 1 for spectral parameters). From left to right: ROSArPSPC spectrum on 1992 May/June (high 
state), XMM-Newton spectrum on 2000 December 26 (high state) and XMM-Newton spectrum on 2001 January 1 (low state). For the XMM-Newton spectra, pn, 
MOSl, and MOS2 data are marked as triangles, soUd squares and circles. 
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Fig. 2.— Long-term lightcurve of the supersoft ULX from 1991 to 2001 (circles: ROSATPSPC; squares: ROSATUBl; triangles: XMM-Newton). The 0.2-2 keV 
luminosity is from spectral fits in Table 1, while the 3a upper limits are assuming a blackbody model with Nh = 10^' cm"^ and kT = 58eV. 
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Fig. 3.— a) Lomb-Scargle periodogram of XMMU J005510.7-373855 as obtained by XMM-Newton on 2001 January. The horizontal dotted Une is the 99.9% 
confidence level; b) Folded Ughtcurve of the "low" state (2001 January) data on a period of 5.4 hr; c) Folded lightcurve of the "high" state (2000 December) data on 
a period of 5.4 hr. Tq of both lightcurves are set at the time of the first data point. 



